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1 A Little History 

Shocks are everywhere in the Universe, and with hindsight it is easy to see 
why. 

In the John Dyson era one might be tempted to ask if indeed there is 
anything else? Such a coherent view of so much of the goings-on between 
the stars has been imparted by John and his colleagues that most interesting 
phenomena now have a Dyson - or Dyson-style - similarity solution to describe 
them. 

But the discovery of shocks in the cold clouds where stars are born caused 
widespread surprise. 

The earliest clear indications came in a rush. As noted in the chapter by 
Tom Ray, it was pointed out bv ISchwart d l|l975|) that the enigmatic Herbig- 
Haro (HH) objects were shock-excited; and in 1976, nearly simultaneously, the 
new techniques of millimetre and infrared astronomy led to important discov- 
eries: 2.3mm CO line p rofiles w ith a width of 150 km s" 1 we r e obse rved by 
iKwan fc Scovillel l)l97fi() and by IZuckerman. Kuiner fc Kuinerl l(l97(tl around 
the position of t he Orion IR cl u ster, a nd in a similar location infrared spec- 
tra obtained by iGautier et alJ l|l976j) showed H2 quadrup ole lines, reveal- 
ing m olecular gas with an apparent temperature of 2000 K IjBeckwith et all 

Ve ry soon a detailed map of shock excited H2 bv lBeckwith. Persson fc Neugebauerl 
lll979h ' appeared, showing two lobes centred around the Becklin Neuge- 
bauer/Kleinmann Low in fra-red cluster; and the important measurement by 
iNadeau fc Geballel l)l979(l of extremely large H2 line- width in the Orion out- 
flow set the scene. 

Al most immediately several i mport ant models were published llHollenbach fc Shulll . 
EES iLonHon. MnOav fc Chul. [l£zl iKwanl ^wanl UajJ demon- 

strated that a hydrodynamic shock travelling faster than 24 km s _1 through 
molecular gas would completely dissociate the H2. This created a major dif- 
ficulty in explaining the wide observed profiles. 
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Th e tw o-fluid magnetic shock f'C type') m odel introduced by iDraind 
(1980) and IDraine. Roberge ~D algarnJ l)l983(l was far less destructive to 
H 2 a nd provided reasonable fits to the line intensities measured at that 
time ijCherneff. Hollenbach fc McKeeL Il983 IDraine fc RobergeL Il98^ . 




Fig. 1. The northern part of the Orion molecular outflow showing the H2 emission 
in a series of bow shocks. The bright tips show up in the light of [Fell], due to strong 
J shocks. [IRIS at AAT picture, courtesy M.G. Burton] 



Meanwhil e the hunt for hot H2 raced on. Many HH objects were observed: 
iMundtJ ( 1987t ) listed fourteen such sources; and the Orion outflow, the bright- 
est source, had been studied in detail. 

A plot of c olumn densities from the near IR lines of H2 measured in the 
Orion outflow ijBrand et all Il988|) demonstrated that a single C type shock 
could not explain the excitation there, and maps of line ratios llBrand et all 
1989) in the outflow, showing little change in excitation with position, further 
clouded the issue. 

As well as these star-forming regions, supernova remnants were mapped in 
shock-excited H2 . The old remnant IC443 shows shocked H2 emission round a 
large part of th e shell, and the line excitation appears strikingly similar to tha t 
found in Orion l|Graham. Wright fc Longmorelll987tlMoorhouse et allll99lh. 
More surprisingly, H2 was found in the heart of the Cvgnus loon ijGraham et all 
1991). 
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In his review iMundd (Il987h suggested one possible explanation for the H 2 
emission often found near, but not exactly at the optical knots in HH objects 
(there was initially some uncertainty because of the low spatial resolution of 
the early IR measurements). Amongst other possibilities he proposed that the 
H2 excitation might come from the flanks of a bow shock, the head of which 
produced the optical emission. 

In this he was following on from the original ISchwartzl l|l978() sugges- 
tion that a bow shock model might be used to explain the range of dif- 
ferent excitations seen in the optical spectra of HH objects. Subsequently, 
lHartmann fc Raymond! l)l984j) proposed that a bow shock model could repli- 
cate the observations of HH1 and HH2; several papers demonst rated that bow 

shock model position/velocity diagram s fitted the observations (Eh oe. Bohm fc Soil 

ll98f)t|Raga fc Bohmlll985tlR.agalll98fii) : and the important paper bv lHartigan. Raymond fc Hartmanr] 
( 1987) explored in detail the excitation and line profile predicted for bow shock 
observations. 

This idea was applied to the observations of the Orion outflow ijSmith fc Brand 
ll990albllSmith. Brand fc MoorhouseLll991b() . demonstrating that single plane 
C-shocks could not explain the observations, but that a C shock propagating 
into material with a very high magnetic field could do so. The magnetic field 
required was extre mely high, implying an u nusual pre-shock magnetic field. 

The review by iDraine fc McKeel l)l993j) summarised the theoretical and 
observational state at the time. 

Before going further, we shall peruse the properties of these shocks. 



2 Shock Basics 

2.1 Hydro dynamic Jumps, Dissociation, Isothermal Shocks, 
Magnetic Fields 

As a prelude, we examine the simplest possible cas e of a steady strong plane 
shock in an ideal gas Ce.g.. lDraine fc McKeel fl993|) . 

Consider the frame in which the shock front is static and the upstream 
gas flows perpendicularly towards it, at a speed V^hock much greater than 
the random thermal motions - hence hypersonically. This is a parallel flow of 
independent particles which encounters a wall of more slowly moving post- 
shock material. Each incoming particle collides elastically (within one or two 
mean free paths) with a post-shock particle, sharing its energy and random- 
izing its direction of travel. Thus, downstream of the shock the gas is hot 
(|fcT ~ 5«iy shock 2 ) and, because of the redirection of velocity, slower. 

It is evident that the shock is a 'sudden' transition in the hydrodynamic 
sense, since hydrodynamics is applicable only to scales much greater than a 
mean free path for elastic collisions; and it remains sudden because to smooth 
itself it would have to propagate pressure waves upstream at the speed of 
sound through gas arriving at a speed much greater than that. 
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At a downstream distance considerably larger than the elastic mean free 
path (in fact of the order of the inelastic mean free path) there has been time 
for collisions to excite the internal states of the gas (which could for example 
lead to ionization and/or dissociation if the shock is strong enough) and to 
change the effective 7 from 5/3 to a smaller number N + 5/N + 3 where N is 
the effective number of internal degrees of freedom of the gas particles. The 
internal states then de-excite by radiating photons and the gas cools. The size 
of this cooling zone depends on the rate of this process, and is large compared 
with the mean free path but usually small compared with the global structure 
of the shocked region. 

More formally, in this frame the mass, momentum and energy are con- 
served through the front. We denote upstream gas by suffix '0', downstream 
by '1'; p is density, p is pressure and v is velocity in this frame, and we as- 
sume that the gas is ideal with constant adiabatic index 7 = 5/3. This value 
corresponds to velocity thermalization by elastic collisions but no internal 
excitation. Hence, the Rankine-Hugoniot shock jump conditions are: 

p\V\ = PoKhock 
Pi + pivf =P0 + Po^shock 2 

7 Pi , 1 2 7 Po . 1 T/ 2 



7-1 Pi 2 1 7-lpo 2 

Further behind the shock when cooling sets in at a net rate C W kg -1 , the 
first two equations, expressing conservation of mass and momentum, continue 
to apply while the third becomes the initial condition for 

V A. { 7 Pl + _ V A = v — = -£ 
dx \7 — 1 pi 2 1 J dx 

in the cooling gas, and 7 may be regarded as a variable to take account of 
internal excitation. 

The various processes giving rise to cooling behind shocks has been 
treated in detail in tw o impo rtant papers bv lHollenbach fc McKeel (flQTflf) and 
iMcKee fc Hollenbachl l|l980ft . Very roughly, in order of decreasing tempera- 
ture, the major coolants are dissociation, collisional excitation of H Lyman 
levels, [01] line emission, H2, H2O (to which any free oxygen is rapidly con- 
verted in the hot gas) and CO line emission. 

Both processes, shocking and cooling, increase the entropy of the gas (the 
first adiabatically) from which we infer the obvious: a shock cannot run in the 
opposite sense. Furthermore, as we saw, a shock travels supersonically with 
respect to the gas ahead. Downstream the gas has to be subsonic with respect 
to the shock so that the high pressure which drives it continues to reach the 
front (the sound speed has been increased to a fraction of V s hock and the exit 
speed has been reduced to a fraction of that speed). 

The quantity being differentiated in the last equation, w, is the stagnation 
enthalpy. If the net cooling function C is expressible as a function of p and T 
then the equation is integrable by straightforward quadrature. 
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Fig. 2. (a) a fast (80 km s _1 ) fully dissociative J shock showing T as a function 
of column density N through the shock. The plateau at 10 4 K is where the atomic 
gas is heated by UV radiation from the front of the shock; the plateau at 4 00 K 
exten ding to N = 10 22 cm -2 is due to H2 reformation ffrom iHollenbach fc McKeeL 
Il989l) . 



If the shock is strong, i.e. the pressure ratio pi/po or equivalently the initial 
Mach number M = V^hock/ (upstream sound speed) is large, then the R-H 
conditions give (with 7 = 1.4, for diatomic molecules without other internal 
states) 

P1/P0 = Vshook/wi = 6, kTx/m = (5/36)V r sh ock 2 , Pi = (5/6)p V s hock 2 - 

As a simple application of these ideas we can show how such shocks have a 
characteristic velocity for the destruction of H2. At the highest temperatures 
the cooling is dominated by dissociation until the temperature drops to a 
value Tii nc at which line cooling dominates and dissociation is much reduced. 
If energy fcldj ss is released per dissociation at a rate R kg" 1 s" 1 then the 
cooling rate is C = kT^ ss R. Each dissociation produces two H atoms with 
number density n(H), so d(n(H) / p) / dt = 2R, whence 

v^-{w + n(H)T diss /2p) = 
ax 

For a strong and initially molecular shock, we neglect the bulk energy in the 
enthalpy, while the pressure is the sum of the partial pressures of atomic 
hydrogen, molecular hydrogen and helium (number fraction y — 0.1), each a 
product of number density and temperature. Define / = n{H)/njj where nu 
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is the number density of H nuclei and is proportional to p; use 7 = 5/3 for 
atoms and 7 = 7/5 for molecules; and note that at the start of the cooling 
zone T = T max and / = 0. Then it is a matter of algebra to show that through 
the cooling flow 

(10^ + 7)(r max -T) 
2Tdiss + 3T 

If we take values of 52, 000 K for Tdi ss and Tii ne = 3500 K, then / becomes 
uni ty - H? is co mpletely dissociated - at T = Tu nc when 14hock=24 km s -1 , 



the lKwanfl|l977|) result. 



In the case where we can treat the cooling zone as small compared to other 
length scales in the region we can lump it together with the shock transition. 
The jump conditions across this 'isothermal shock' (for simplicity it is assumed 
that initial and final temperatures - and internal states - of the gas are the 
same: hence the name) are 

Pi/ Pa = M? so , pi = po^shock 2 - 

Cooling allows the density to reach very high values set by the 'isothermal 
Mach number' M; so = 14hock/ \/kTo/m. In the cooling zone the pressure has 
increased by only 17%. 

If a magnetic field is flux-frozen to the gas (i.e. a fully ionized plasma) and 
the field is directed at right angles to the shock velocity, then pressure is the 
sum of that of plasma and magnetic field. The magnetic pressure is B 2 /2/xq- 
The effective mass density of the magnetic field for non-rclativistic shocks is 
negligible, and the shock is assumed to be strong. The momentum equation 
is then 

p + B 2 /2p Q + pv 2 = B 2 /2p + PoV shock 2 

while mass conservation becomes pv = po^hock- Since flux freezing is assumed, 
B oc p oc v^ 1 . Denote the Alfven speed by Va (V| = B 2 /2p p), and the 
Alfvenic Mach number by Ma = V^hock/VA.o- 

Writing v = v/V s hock and dividing through by poKihock 2 the momentum 
equation becomes 

(p/p V shock 2 ) + (2MIV 2 )- 1 +v = {2M 2 A )- X + 1 

and if downstream the gas is cool enough that magnetic pr essure dominates 
(for more general conditions see iRoberee k Drainel ll99C|) then the down- 
stream value of u/V^hock is 



(1 + y / l + 8Ml)/4M 



For magnetic pressure domination downstream (i.e. the final Alfven speed 
greater than the sound speed, assumed equal to the initial sound speed for 
simplicity), assuming as we have that the Mach number M is large, 



M A < V2M 2 . 
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In such shocks the cooling zone is supported by the magnetic field stress at 
a density not hugely different from the initial density (thus slowing cooling 
processes in comparison with that in the high density achieved in the absence 
of magnetic field). 

The detailed structure of J-shocks is modified by the possible changes 
of state of the material, and its effect on the surroundings. If the shock is 
strong enough, radiation generated in this region will propagate ahead of the 
shock and may alter its state depending on shock velocity. And since hydrogen 
molecule re-formation is a slow process, dissociated molecules will not reform 
in the close vicinity of the sh ock, and a separate re-formation zone may occur 
llHollenbach fc McKeeL llQSfll) . 

2.2 C Shocks 

A quarter of a century ago, following an idea bv lMullanHl97lft . lDrainel Jl98f]l) 
pointed out that, since molecular clouds are slightly ionized and are permeated 
by magnetic fields, a different kind of shock is possible and might be prevalent. 
This he named the C shock, using the name J shock for the hydrodynamic 
discontinuity just described. 

The C shock can be envisaged in a two fluid continuum. One of the fluids is 
the molecular neutral gas representing most of the mass, while the other fluid 
is a conducting magnetized plasma consisting of the tiny fraction (typically 
10~ 7 — 10~ 5 ) of particles ionized by cosmic rays plus the magnetic field in- 
ferred from observation, in which the plasma is presumed to be flux-frozen. If 
a smooth pressure gradient is generated perpendicular to the magnetic field it 
will drive a compressive magnetohydrodynamic wave through the plasma. This 
wave, if unencumbered by interaction with neutral particles, would rapidly 
steepen into a shock travelling faster than the magnetosonic speed (in these 
cases roughly equal to the ionic Alfven speed, much greater than either the 
thermal sound speed in the neutral gas or the Alfven speed Ma of the whole 
medium treated as a single magnetized fluid). But by virtue of gas viscosity 
the magnetic pressure is transferred slowly to acceleration of the neutral ma- 
terial until far downstream the two components have the same velocity. The 
viscous drag can ensure that the plasma pressure wave will not accelerate to 
a shock state. For a wide range of conditions (magnetic field strength, state of 
ionization) the neutral gas is accelerated slowly, and remains supersonic with 
respect to the structure throughout, without passing through a hydrodynamic 
shock. The whole process reaches a stationary state for suitable initial condi- 
tions far upstream and boundary conditions downstream. The relatively wide 
region where this transition occurs is the C shock. In this region the neutral 
gas is heated by the viscous friction, and this region radiates. 

The structure of this type of shock is most easily considered using the 
mass and momentum conservation equations introduced in the discussion of 
J shocks, applied separately to the neutrals and to the magnetized ion fluid. 
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C shock 




Fig. 3. Sketch, in the shock frame, of a J shock showing T and v, the dashed 
line corresponding to the effect of the presence of a magnetic field, and a C shock, 
showing the slab of gas at a relatively uniform temperature. 



This time the neutral fluid has an extra force on it caused by friction 
with the ions, and vice versa. But in the summed momentum equation these 
terms cancel, and we can use the equation neglecting these forces. Generalizing 
the result for a magnetized one-fluid shock above by noting that flux freezing 
applies to the ions only, and denoting ion properties by subscript i and neutrals 
by n, 

(Pn/pn,oKhock 2 ) + {2M\v1)- 1 +V n = (2MI)- 1 + 1 

whence we find a relationship between the ion and neutral velocities. 

If the cooling is such that thermal pressure is unimportant, then approxi- 
mately 

v n = l- {vf - 1)/2M|. 

Far downstream where v$ = v n we have the same conditions as the cold 
magnetized shock discussed above. 

The slip speed, that is Vi — v n , determines the amount of frictional force 
and fri ctional heating in the flow. In the sim ple case here we find the maximum 
value llSmitli. Brand I fc Moorhou ^el ll 991 hh 



= l + l/2Ml-3/2M 2 A 



Consider the length scales in this shock structure. The frictional body 
force (per unit volume) imposed on the neutrals by the ions, and vice versa, 
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Fig. 4. (a) A low Ma C shock showing Vi and v n changing in parallel, (b) A high Ma 
C sho ck sowing the rapid change in Vj an d the consequent large velocity difference 
(from ISmith. Brand fc MoorhouseL Il991al^ . 



has magnitude F = p n Pi{w){v n — Ui)/(m, + m n ) ijDraine fc McKeei . Il993(l . 
where (crv) is the ion-neutral elastic collision rate, and subscripts n, i have the 
obvious meaning. The equations of motion for ion fluid (neglecting ion thermal 
and ram pressure) and neutrals (neglecting thermal pressure) respectively are 

-(-)=» 

dz \2fi a J 

If we define Lq — 14hock/"-i.o(o'w), roughly the initial mean free path of neutrals 
through the ions times Mach number, and define also the ion length Li = 
Vi I {dvi / ' dz) and neutral length L n = v n /(dv n /dz), the first equation above 
gives Li I Lq — M^ 2 and we also find from the (vi,v n ) relation that L n /Li = 
M A . 

These relations for slip speed and length scales show the following results. 
If Ma is small the slip speed is a small fraction of the shock velocity. Such soft 
shocks can accelerate the neutrals without getting them very hot. The change 
in the ion velocity closely parallels that in the neutrals. The length of the slip 
zone is proportional to 14hock, and so can accelerate the radiating molecules 
to quite large velocities. 

On the other hand a high value of Ma will create a flow in which the ions 
rapidly approach their final speed before the neutrals have been significantly 
decelerated, exposing the neutrals to a high rate of high velocity collisions. 
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There are two major consequences. The heating rate is increased, and the 
gas may become sufficiently hot that it becomes subsonic - which engenders 
a hydrodynamic jump (the neglected neutral thermal pressure term in the 
equation above becomes important). And in these conditions at a shock ve- 
locity of 40-50 km s _1 the neutrals can be ionized by single collisions, leading 
to an ionization run-away as the increasing ion fraction reduces Lo- 
in any case if the length of the drag zone is too short compared with the 
cooling length, the sound speed will increase and the gas may become subsonic 
through a hydrodynamic jump. Thus, for a limited range of conditions a C 
shock structure will contain a J shock which may provide much of the cooling 
via high excitation lines but still produce C shock-like emission of H2 for 
example. 

The heated region behind each kind of shock is shown schematically in 
Fig.H 

It is evident that, in contrast to J shocks, C shocks create a zone that very 
roughly corresponds to a slab of gas at a given temperature, and that therefore 
the emission will reflect that status. That this is not shown in m any of the 
observations is the reason for proposing bow sho cks as emitters l)Schwartzl 
Il975l lHartmann & Raymond Il984t iMundtl fl987^ . 

2.3 Bow Shock Structure 

There are few instances of observations of planar shocks. Usually the next best 
guess is that one is observing a shock driven into a bowed shape by a super- 
sonic gas jet or supcrsonically moving dense clump. The first models of Herbig 
Haro objects were of this nature, and expose the essential complications. 

The range of excitations and breadth of velocity profiles seen in quite small 
projected areas on the sky led to the idea that one is often observing some 
or indeed all of a bow shock structure. In some special cases the bow is quite 
discernible. 

A simple but perhaps suspect way in which to create models of such 
regions is to use oblique planar models, and patchwork these onto the 
surface of a pre-decided bow geometry . This technique was pioneered by 
lHartigan. Raymond &: Hartmarml Jl987l) for opt ically emitting bows, and ap- 
plied to H 2 llflmith. Brand fc Moorhon sel ll 991 1 ah to explain the extreme con- 
ditions at the Peak 1 in the Orion outflow. 

The strongest section of the shock is clearly at the head of the bow. Here, 
the shock may be strong enough to ionize and dissociate the molecules, and is 
certainly a J-shock. Further round the side of the bow the shock may still be 
a J-shock but not completely dissociating the molecul es, while further round, 
the shock becomes a C-shock ijSmith fc P,randUl990rh . 

Clearly, even if the initial magnetic field had a simple geometry, the varying 
angle of incidence at each part of the shock makes such a structure hard to 
model. The extended length scales and subtleties of structure of oblique C 
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shocks l|Warclle fc Drainel Il987|) make this region the least straightforward 
part of the modellin g process. 

Notwithstanding. lSmith. Brand fc Moorhousd l|l991aj) were able to repro- 
duce both the excitation (up to very high energies) and the extreme breadth 
of the lines with a very low Alfven Mach number and very fast shock. The lin- 
gering problem with this explanation has always been seen to be the extremely 
high pre-shock magnetic field required. 

The technique is to specify an axially symmetric bow shape R = f(z) (z 
is distance along the axis; R is perpendicular distance from axis of the bow at 
z) supposedly created by a uniform supersonic flow, velocity Kjhock, hitting a 
small obstacle, or equivalently by such an obstacle - or indeed the end of a 
jet - ramming supersonically into a uniform medium. If dz/dR = tan^ the 
effective shock speed at a point in the bow is V s hock cos -0 and the shocked 
material has, in addition to its post-shock behaviour perpendicular to the 
shock, a component of velocity K,hock sin ip tangential to the bow. A line of 
sight and a projected area around it are assumed to model an observation. At 
each point of the bow within that area a line intensity (per unit area of shock) 
through the shock flow is calculated, as are the projection factor (secant of 
angle between line of sight and shock normal), and the velocity component 
along line of sight. From this the line profile is found, the shape and intensity 
of the line being diagnostic of the shock properties. 

3 Models and Observations 
3.1 Issues 

Much of the theoretical work until recently has been based on the assumption 
that the flow is stationary, and that timescales for chemical evolution were 
long, and timescales for population of internal states were short, and that 
upstream conditions are static. Even now the modelling of shocks, including 
chemistry, grain physics and electrodynamics and the effects of the dust on 
the chemistry and dynamics, and obliquity progresses inevitably slowly. The 
fact that post-shock turbulence will occur in many of the flow structures that 
have been measured has led to some pioneering studies of the process, but 
much remains to be done. 

The question of how to extend a two-fluid approximation to account for 
the behaviour and effects of the dust grain population has been addressed 
several times, and because of the still uncertain properties of grains is still an 
open question. 

Even in plane steady normal shocks several instabilities have been discov- 
ered. Are there more? Is it possible to make useful predictions from unstable 
shocks? 

On the observational side, the ISO measurements, opening wider the near 
to far IR spectral region, have produced a plethora of new results. But there 
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is still the usual urgent need of better spatial resolution (and of course signal- 
to-noise ratio) to compensate for our inability to see round the side of these 
complex regions. 

Grain Charge and Drag 

The picture so far of the C-shock is of a two- fluid process, but the presence 
of charged dust grai ns can s i gnific antly alter the picture, particularly in cold 
high density regions l|DraineL Il98(l) . Severa l papers investigate the tricky issue 
of grains as a subspecies in the shock flow llPilipp Hartauist fc Havnest ll99Ct 
Caselli. Hartauist fc Havneslll997tlWardleUl998llFlower fc Pineau des Foretst 
20031) . Also see the chapter by Tom Hartquist and Ove Havnes. 

The first issue is that the grains may be charged, and can be the dominant 
ionized species pr oviding the drag on neutrals . This will occur preferentially 
at high densities. Pili pp. Hartauist fc Havnesl 1)199(1) found that the dust in 
high density shocks was constrained by charge separation to move at a speed 
intermediate between that of th e ions and of th e neutrals, with the net effect of 
steepening the shock structure. IWardlel l|1998h examined the possible flows in 
great detail, and pointed out that the effect is accentuated by taking the out- 
of- plane deviations of the magnetic field into account, thus finding even steeper 
(and by implication hotter) shock structures. These effects are greatest at high 
density and in the presence of sm all grains (which for a given grain mass will 
give greater dragV TWardkj l|l998|) noted that, since most C shocks will give 
temperatures between 1000K and 2000K, and molecules at that temperature 
will radiate most of the shock energy, there may be no dominant observed 
effect in changing the detailed structure of the C shock. Conversely, he pointed 
out that reliable diagnostics will rely on getting the physical processes correct 
in detail. 

iFlower fc Pineau des Foretsl l)2003h addressed some of the physical and 
chemical detail. They show that large grains rapidly charge by electron at- 
tachment in the shock flow and become attached to the magnetic field, and 
that the critical velocity ( at which H dissociation leads to breakdown to a J 
shock) i s sufficiently high llLe Bourlot et "ail . 12002^1 that refractory grains can 
shatter. IChapman fc Wardlel l|2006) examined the effects of obliquity of the 
magnetic field, and emphasi zed the significance of the grain physi cs in de- 
termining the flow structure. ICiolek. Roberge. fc Mouschoviasl (|2004f) pointed 
out that the degree to which the grains attach to the magnetic field deter- 
mines the ionic magnetoionic sound speed, and if a significant fraction of the 
total grain mass is thus attached, the ionic material will shock, destroying the 
C shock. 

Grain Destruction 

In all of the calculations referred to above there is the issue of grain destruction 
(or more positively element release) cither by sputtering or shattering of dust 
grains. 
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Calculations by iDraind Jl995l) and iFlower fc Pineau des Foretsl l|l995|) 

have shown that C shocks can de stroy dust grain mantles and cores, and 
Caselli. Hartauist fc Havnesl l|l99lfl have shown that such shocks may release 
an abundance of Si consistent with observations o f shocked regions, and many 
orders of magnitude above the average abundance. ICaselli. Hartauist fc Havnes 
ljl997h emphasized that destruction is more efficient in oblique shocks than 
in perpendicular shocks. This process may be the dominant producer of SiO 
seen in high density star fo rming regions and indicated in the observations of 
Jimenez- Serra et alJ ll2005l). for examp le. 

" lMoweT^TF^hie^mdesForetsI l)2003|) demonstrated that a specific C shock 
travelling at 50 km s _1 would shatter over half of the amorphous carbon grains 
present in the model. 

All of these findings together imply that the behaviour of C shocks is 
dependent on dealing with chemistry, dust grain physics and MHD phenom- 
ena in great detail (since some of the necessary grain physics remains to be 
determined, this is a worry). 



Chemistry and Time Dependence 

Another issue which has been clarified by the refinement of shock models is 
that the processes determining the shocks have timescales which are com- 
mensurate with shock passage times, but also may be comparable with the 
time since the shock was created. This raises the possibility for the pessimist 
that every shock observed is unique (or rather that the number of determin- 
ing parameters, now including initial conditions and time since creation, is 
far greater than observations can encompass). At any rate it emphasizes how 
important the detail can be in determining the overall effect. 

Further, if the pressure driving the shock changes, or the medium into 
which it propagates has large gradients, or is lumpy, externally imposed time 
dependent effects come into play. 

An example of what time-dependent driving can do is shown bv lLim et al.l 
They demonstrated that a shocked layer driven by a slowly accelerat- 
ing piston can accumulate molecular gas without destroying the H2, producing 
column densities and velocities of gas comparable with those observed in the 
Orion outflow in an acceleration period of tens of years, and protecting the 
molecular layer from the hot gas at the shock by an intermediate layer of H 
at oms. 

iGirart et alJ <)2005[) have demonstrated that near HH2, and caused by it, 
a wide variety of different excitation phenomena can be created by a shock in 
a non-uniform medium. 

IChieze. Pineau des Forets fc Flowerl ( 1998^ demonstrated that the timescale 



for an initial disturbance in a dark cloud to reach steady state is so long that 
it may exceed the lifetime of the outflow causing the shock. They showed fur- 
ther that during the evolution to a C shock there will be a J-shock embedded, 
and that the combined radiation from this region and from the rest of the 
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(C type) flow may explain the H 2 observations in some sources where several 
different excitation temperatures are measured. They emphasized again the 
necessity of calculating in detail the degree of ionization in particular. 

These more detailed studi es also reveal t hat e ven if a J shock occurs a ll 
H 2 emission may not be lost! iFlower et all l)2003l) and ICabrit et all l)2004j) l 
demonstrated that for fully dissociating shocks of moderate speed an H2 emis- 
sion spectrum can be produced which can account for some of the observations. 

The issue of time dependence is intimately tied up with instabilities. The 
most marked(!) instability is the lWardlel l)l99d) instability of plane C shocks, 
in which, rather like the Parker instability, ripples on the lateral magnetic 
field lines accumulate a higher density of ions and are dragged further down- 
stream, trapping higher density gas which cools rapidly. The as tonishing thin g 
about this instability, which also occurs in oblique C shocks l|Wardlelll99l . 
is that it has little effect on the overall output fr om the shock, for reasons 
touched on before: the hot gas does the radiating ( Mac Low fc SmithL Il997t 
iNeufeld fc StoneUl997lk 

Many of the developing models reveal one dimensional chaotic or bounc- 
ing instabilities, due for example t o shortened cooling times in disso c iation 
or ionization zones behind shocks l(Lim et al 1 l2002l lSrnith fc Rosen! bon,^ 
iLesaffre et aTll2004|) . "What will happen in three dimensional models?" is an 
open question. 



Turbulence and Turbulent Mixing and Other Effects 

I Canto fc Ragall|l99l[K in a gr ound-breaking paper, in vestigated heating effects 
of a turbulent mixing layer. iTavlor fc Raeal (1995) showed that such zones 
could mimic post-shock flow s in som e res pects. Amongs t other s, IPavlovski et al. 
l|2002l) . lElmegreen fc Scald l(2004) . and iHeitsch et alJ l|2005h have begun to 
model more general cloud turbulence as a source of structure and radiation. 
This work is still at an exploratory stage, and in future may converge with 
that on shock codes. 



3.2 Current Models and Observations 

In the interstellar medium in our Galaxy we expect that the cooling zone 
in most observed shocks will be ~ 10~ 5 pc in the case of J-shocks and ~ 
10 -2 pc for C-shocks in very round numbers. This means that it is unlikely 
that J-shocks will be resolved, while C-shocks can be resolved with current 
equipment. 

T he arrival of la rge masses of ISO satellite data has enlivened this 
area l|DishoeckL Eool . 



Modelling is rapidly b ecoming more sophisticated, wit h attention being 
paid to optimum methods IjFallet l2003i iLesaffre et all l2004j) . and more of the 
required physics (g rains, chemistry, time dependence) is being included. What 
has become clear jChieze. Pineau des Forets fc Flower! Il998l IFlower et all 
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l2003HLesaffre et all Eool ICabrit et~afl 12004^ is that rugged predictions re- 
quire very complete modelling. 

Now a few examples are considered. 

The extensi ve line list of H? from HH43, observed bv lGiannini et all {2002) , 
has been fit by iFlower et all l|2003j) by either a 25 km s" 1 J sho c k or a n 80 
km s -1 C shock! The same mode llers fit data by IWright et alJ l|l996h and 
iFroebrich. Smith fc Eisloffell j2002j) of Ceph A West with a model of a J shock 
wit h a C precursor t ypical of early evolution. 

lO'Connell et alJ (2005) modelled the H2 emission structures seen in the 
HH2 11 outflow using C bow shocks. 

ISnell et all feOOfiji examined the supernova remnant IC443 using the Sub- 
millimeter Wave Astronomy Satellite (SWAS) to observe H2O lines. They 
concluded that the data are consistent with a combination of fast J shocks 
and slow (J or C) shocks, perhaps likely in a clumpy medium. 

The Orion outflow, pictured in part in Fig. ^ has been a test-bed and 
co nundrum throughout the entire development of this area. 

iLe Bourlot, et alJ ( 20021) produced the H? excitation diagra m which matches 



the observed data of iRosenthal. Bertoldi fc Draoatd l|200f]|) by superposing 
two C shocks with velocities as high as 40 k m s" 1 and 60 km s" 1 . 

The discovery of the 'bullets ' or fingers lAllen fc Burtonl dl993l) fsee also 

McCaughrcan & Mac LowMl996l) suggested an explosive event. and lStone. Xu fc Mundvl 
( 1995) provided an explanation for the major structural features (in particular 
the fingers) by having a time variable outflow, from IRc2 or the BN object, 
which sweeps up a shell and then accelerates it whereupon it fragments via 
Rayleigh- Taylor instability. 

There has been a great deal of progress in establishing the values for 
important rates in the microphysics and chemistry, in developing models that 
can cope with the subtleties of grain dynamics and time-dependent chemistry, 
and in the observations (particularly thanks to ISO in the mid IR); but an 
enormous amount remains to be done. 
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